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Abstract 

We study partitions of the two-dimensional flat torus (R/Z) x (R/5Z) into k domains, with b a 
real parameter in (0,1] and k an integer. We look for partitions which minimize the energy, defined 
as the largest first eigenvalue of the Dirichlet Laplacian on the domains of the partition. We are in 
particular interested in the way these minimal partitions change when b is varied. We present here 
an improvement, when k is odd, of the results on transition values of b established by B. Helffer and 
T. Hoffmann-Ostenhof (2014) in [14] and state a conjecture on those transition values. We establish 
an improved upper bound of the minimal energy by explicitly constructing hexagonal tilings of the 
torus. These tilings are close to the partitions obtained from a systematic numerical study based 
on an optimization algorithm adapted from B. Bourdin, D. Bucur, and E. Oudet (2009) in [4]. 
These numerical results also support our conjecture concerning the transition values and give better 
estimates near those transition values. 

MSC classification. Primary 49Q10; Secondary 35J05, 65K10, 65N06, 65N25. 

Keywords. Minimal partitions, shape optimization, Dirichlet-Laplacian eigenvalues, finite differ¬ 
ence method, projected gradient algorithm. 


1 Introduction 

1.1 Motivation 

Optimal partition problems are a field of shape optimization which has recently generated much interest, 
see e.g. laElIHlIls]. They are connected to nodal sets of Laplacian eigenfunctions (see e.g. m) and to 
steady states for competition-diffusion systems of partial differential equations (see e.g. M)- This paper 
focus on a particular type of problem, studied by B. Helffer, T. Hoffmann-Ostenhof, and S. Terracini in 

m- 

Let us describe the general setting. In the following, M is a compact, two-dimensional, Riemannian 
manifold without boundary, and k an integer, k > 2 . All along this paper, we consider /^-partitions of 
M in the following sense. 

Definition 1.1. We call /^-partition (or simply partition^ a finite family V = of open, 

connected, and mutually disjoint subsets of M (called domains of the partition). This partition is strong 
if M = IJ2=i case, we can define the boundary of V as N(fD) = |Ji=i then say 

that V is regular if N(fD) is locally a regular curve, except at a finite number of singular points, where 
a finite number of half-curves meet with equal angles (the ’equal angle meeting property’). 

We denote by the set of all k-partitions. 

For any /c-partition V G , we define its energy by 

Ak{V) = max Ai(A), (1-1) 

l<i<k 
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where, for any open set cj C M, we denote by (Aj(cj))j>i the eigenvalues of the Laplacian on cj, with 
Dirichlet boundary conditions, arranged in non-decreasing order and counted with multiplicities. The 
optimization problem we consider here is to study 

2k{M)= inf Afe(P). (1.2) 

We say that a partition P* G is minimal when Kk{V'') = £/c(M). The existence and regularity of 
minimal partitions have been established in [71 HOl El [15] , and we have 

Theorem 1.2. Let V = {Di)i<i<k be a minimal k-partition. Up to zero eapaeity sets, V is strong and 
regular. Furthermore, V is equispectral, that is to say A/c(T>) = Ai(T)j) for any I < j < k . 

In the subsequent paper mi, B. Helffer and T. Hoffmann-Ostenhof considered the case where M is 
the two-dimensional flat torus T(l, 6) = (M/Z) x (M/6Z), with b a parameter in (0,1] . We also consider 
this situation in the present paper. In particular, we are concerned with finding numerically the minimal 
partitions of T(l,6), for a range of values of b. Minimal partition problems on surfaces have already 
been investigated numerically by several authors. The paper [4] treats the case of the torus T(l,l), 
and was the inspiration for our numerical study, while m considers several surfaces embedded in , 
including the sphere. A more complete study of the sphere, using a different method, is presented in [1]. 
In all these papers, the energy to be minimized is the sum of the eigenvalues rather than the maximum. 


1.2 Nodal partitions 

We will exploit the connection, shown in m, between minimal partitions and nodal domains of eigen¬ 
functions of the Laplacian. Let us recall some relevant definitions and results. With an eigenfunction u 
of the Laplacian on M , we associate the nodal domains which are the connected components, denoted by 
Di , of M \ N{u) with N{u) = {x G M; u{x) = 0} . The number jii{u) of connected components is finite. 
According to classical results on the regularity of the nodal set, the family of the nodal 

domains is a strong regular /i(i4)-partition in the sense of Definition |1.1[ We call it the nodal partition of 
u . Let us note that, for a nodal partition, an even number of half-curves meet at each singular point in 
the boundary. This is in contrast with the situation for minimal partitions, where there is no constraint 
on this number. 

A famous result, proved by R. Courant m, states that if u is an eigenfunction associated with Xj{M) , 
u{u) < j • Following [15], we introduce a new definition for the case of equality. 

Definition 1.3. An eigenfunetion u of the Laplaeian, assoeiated with the eigenvalue X, is said to be 
Courant-sharp if ii{u) = min{^; A^(M) = A} . 


We can now state a result of m that links minimal and nodal partitions. 


Theorem 1.4. We have Xk{M) < £/c(M), and therefore the nodal partition of a Courant-sharp eigen¬ 
funetion is minimal. Conversely, if the nodal partition of some eigenfunetion is minimal, this eigenfune¬ 
tion is Courant-sharp. Finally, if 2k{M) = Xk{M), all minimal k-partitions are nodal. 


Remark 1.5. Theorem l.f implies that 2k{M) 
and 2k{M) > Xk{M) otherwise. 


Xk{M) if there exists a nodal minimal k-partition. 


Let us note that eigenfunctions associated with A 2 (M) always have two nodal domains. Therefore, 
according to Theorem |1.4[ the notion of minimal 2-partition coincides with the notion of nodal partition 
of a second eigenfunction. The problem of finding minimal /c-partitions becomes interesting for /c > 3, 
which we will assume in the rest of the paper. 


1.3 Summary of the results 

To use the nodal partitions, let us give explicitly the eigenvalues of the Laplacian on the general torus 

T(a, b) = (M/aZ) x (IR/6Z) , with 0 < 6 < a. 

The Laplacian on T(a, b) is unitarily equivalent to the Laplacian on the rectangle R(a, b) = (0, a) x (0,6) 
with periodic boundary conditions, and its spectrum can be computed by separation of variables. 
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Proposition 1.6. The eigenvalues of the Laplaeian on T(a,6) are 


/^2 ^ 2 \ 

Am,n{a, &) =47r^ (-^ +-p-j , with m, n e No . 

Let us recall that, according to m, the only non-constant Conrant- shar p eigenfunctions for the torus 
T(l, 1) are associated with A2(T(1,1)) = 47r^ . According to Theorem 1.4, this implies that, as soon as 


k>3, a minimal /^-partition of T(l, 1) is not nodal and we have to find new candidates. 
Let us introduce a particular partition of T(a, b) into vertical strips. 

Definition 1.7. We denote by Vk{a^h) the k-partition o/T(a, 6) with domains 

'i-1 


A = 


k 


-a, -a ] X (0,6) , 


fori = l,...,k. 


We have Kk{Vk{a^b)) = k'^'K^jo? ^ and any partition obtained from Vk{a^b) by a translation has 
the same energy. Let us note that when k is even, Vk{a^b) is the nodal partition of the eigenfunction 
{x^y) ^ sin(/c7rx/a), associated with the eigenvalue A/j./ 2 ,o (<^5 ^)- 

Let us now focus on the torus T(l,6), with 6 G (0,1] . Following [14], we want to know for which 
values of b the partition P/e(l, 6) is minimal. Then we define the transition value b^ by 


bk = sup{6 > 0 ; T>/c(l, b) is a minimal /c-partition of T(l, 6)} . 


(1.3) 


This notion is well adapted since T>/c(l, b) is a minimal /^-partition of T(l, 6), for any b G (0, 6/^], as it will 


be established in Proposition 2.1 


When k is even, a direct application of Theorem 1.4 gives us the transition value (this result and its 
detailed proof can be found in El)- 


Proposition 1.8. If k is even, then b^ = 2/k . Furthermore, ifb < 2/k, 'Dj^{l,b) is, up to a translation, 
the only minimal k-partition o/T(l,6) . 


Let us note that if k is even and 6 = 6/c, P/c(l, 6) is no longer the only minimal partition of T(l, 6) 
up to a translation, due to the multiplicity of the eigenvalue A/c(T(l, 6)). We will see this in more detail 
in Section [2~3l 


When k is odd, m Theorem 1.1] proves that bk > 1/k. Before stating our improvement of this 
estimate, let us introduce some notation. For b G (0,1] , we consider the infinite strip = M x (0, b) 
and we define 

bl = sup {b G (0,1] ; j{b) > , with j{b) = inf Ai(fl). (1.4) 

QcSb,\Q\<b 

Theorem 1.9. If k is odd, then b^ > b^ > 1/k. 

We are also interested in obtaining upper bounds of £/c(T(l,6)) for 6 G (0,1] . Since Ak{Vk{l,b)) = 
, we always have £/c(T(l,6)) < . For some values of the parameters, we construct hexagonal 

partitions which give us an improved upper bound. 

Theorem 1.10. For /c G {3,4, 5} , there exists b^ G (0,1) sueh that, for any b G (6^, 1] , there exists 
a tiling o/T(l,6) by k hexagons that satisfies the equal angle meeting property. We denote by H/c(6) the 
eorresponding tiling domain, and we have 


£fe{T(l, b)) < min Ai(Hfe(6))), V6 e (6^, 1]. 

More explieitly, we ean ehoose 

^3 = ^ — 0.396, 64 = 0.289 < 64 = ^ , and 65 = ^ 0.233. 


2^3 


36 
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The paper is organized as follows. In Section we analyse the transition values and give a proof 
of Theorem |1.9[ We also set out some arguments to conjecture the transition values and for such 6, 
we present some candidates to be minimal partitions, which are obtained from eigenfunctions of the 
Laplacian on a covering of the torus. In Section we describe our numerical method, which is based on 
the work of B. Bourdin, D. Bucur, and E. Oudet (2009) in [4]. For a fixed /c G {3,4, 5}, we compute 
candidates to be a minimal /c-partition of T(l, 6). These computations suggest the existence of hexagonal 
tilings of a specific type. In Section we construct explicitly these tilings and compute their energy. 
This improves the previously known bounds of the minimal energy. Near the conjectured transition 
values and near 6 = 1, the numerical simulations give us better candidates. 


2 Transitions between different types of minimal partitions 

2.1 Transition value 

Let us first show that the notion of transition value introduced in ( |1.3| ) is well behaved. 

Proposition 2.1. For any b G (0,6/^], Vf^{l^b) is a minimal k-partition ofT{l,b) . 

The proof is a direct consequence of the following properties of 6 £/c(T(l, 6)): 

Proposition 2.2. The function 6 £/c(T(l, 6)), defined on (0,1], is continuous and non-increasing. 

Proof Let us pick b and b' in (0,1], with b' < b. We define a mapping F from T(l, b') to T(l, b) by 

F: T(l,6') ^ T(l,6) 

{x,y) {x,^y) . 

If uj is an open set in T(l,6'), F{uj) is an open set in T(l,6). By a direct estimate of the Rayleigh 
quotients, using the change of variable defined by F, we obtain 

Ai(FH)<AiH<^Ai(FH). (2.1) 

Let us now consider a minimal partition V' = (F-)i<^</c of T(l,6'). We define the partition F{fD') = 
{F{D[))i<i<k of T(l, b) . According to the inequality on the left in ( |2.1[ ), we have 

b)) < Ak{F{V')) < Ak{V') = £,(T(1,6')). 


Since b and b' are arbitrary, this establishes monotonicity. Let us on the other hand consider a minimal 
partition V = (A)i<i</c of T(l,6), and define the partition F~^{'D) = (F“^(F^))i<^</c of T(l,6'). 
According to the inequality on the right in (2.1), we have 

i:fe(T(l, 6')) < MF-\V)) < ^Ak{V) = b)). 


Since b and b' are arbitrary, this establishes continuity. 


□ 


2.2 Proof of Theorem 11.91 

We break down the proof into several lemmas. The structure of the argument follows closely la, with 
two main changes. First, we have imposed the constraint |ci;| < 6 in the auxiliary optimization problem 
( |1.4[ ), in addition to the inclusion constraint u; C . This allows us to exclude the existence of domains 
homeomorphic to a disk for a larger range of values of 6, as discussed in Lemmas |2.3| and |2.4[ We 


have also used the pair-symmetric structure of the lifted partition to obtain a better lower bound of the 


energy, as seen in Lemmas 2.5 and|2.6 


The definition of the optimization problem (1.4) implies that b t-G j{b) is non-increasing with respect 
to b . According to the Faber-Krahn inequality, the disk is a minimizer for sufficiently large b . For a 
small b , the following lower bound, deduced from the Poincare inequality on , gives more information 
(see [H]): 

TT^ 

J(‘) a ^ ^ 
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It can be shown, using a concentration-compactness result for shapes proved by D. Bucur (see 0 ), that 
there exists a quasi-open minimizer for Problem (1.4), which satisfies |f7*| = b. This allows us to 
obtain an estimate of . 


Lemma 2.3. We have 


l<bl< 


1 


Proof. An elementary argument using monotonicity, similar to the proof of Proposition | 2 . 2 [ shows that 
the function b ^ j(b) is continuous. Furthermore, in the case b = l//c, we have seen that there exists a 
quasi-open set fP C S 5 such that Xi{PP) = and |f]*| = 1/k. This condition on the area implies 

in particular that , by strict monotonicity of the first Dirichlet-Laplacian 

eigenvalue. We conclude that > 1/k. 

To obtain the upper bound, let us consider the rectangle 


R(l, 6 ) = ( 0 ,l)x( 0 , 6 ) . 

We have R(l, 6 ) C S 5 and |R(1,6)| =b. If b = — 1, we have Ai(R(l, 6 )) = . Furthermore, we 

know that R(l, 6 ) is not minimal, since the normal derivative of a first eigenfunction is not constant on 
the free boundary. We conclude that 6 ^ < 1/ \/k‘^ — 1. □ 

We now prove a topological property of minimal partitions. 

Lemma 2.4. Ifb<b\, a minimal partition o/T(l, 6 ) has no domain homeomorphie to a disk. 

Proof. Let us assume that some minimal partition V has one such domain, that we denote by D. We 
consider the following covering of T(l, 6 ) by the plane : 

Boo: ^ T(l, 6 ) 

(x^y) 1 -^ {x mod 1 ,^ mod b). 

Let Dq be one of the connected components of (D). It is homeomorphie to a disk, and \Do \ = \D\ < b. 
Furthermore, since T(l, 6 ) is of width b , for any xq G M, the total length of the vertical slice at xq , that 
is to say of the set {y ; (xq, y) e Dq} C R ^ is smaller than b. Let us call Dq the Steiner symmetrization 
of Do with respect to the line {{x^y); y = 1/2} . It has the same area as Dq , and it is contained in , 
according to the geometrical property mentioned above. Since Steiner symmetrization does not increase 
the first eigenvalue, we obtain 

j{b) < Ai(D^) < Ai(Do) = Ai(Z)) = AkiV). 

Since 6 < , j{b) > = A/. b )), contradicting the minimality of V . □ 

Following we consider the torus 

T(2,26) = (M/2Z) X (M/26Z) 

equipped with the natural projection map 114 • {x^y) ^ {x mod l^y mod b) from T(2,26) to T(l, 6 ). 
Since every point of T(l, b) has four antecedents by 114 , T( 2 , 2b) can be seen as a four-sheeted covering 
of T(l, b). The pull-back Ilf^{D) of a connected open set D in T(l, b) is an open set in T(2, 2b) having 
at most four connected components. With any /^-partition V = (D^)i<^</c of T(l, 6 ), we can therefore 
associate the partition whose domains are all the connected components of all the sets llf^{Di ), for 
i G {1, ..., /c} . We call it the partition lifted from V and denote it by Ilf^{V). It is a regular l- 
partition of T(2, 2b ), with k < £ < Ak, and it has the same energy as V. 

Let us now define the following mapping on T( 2 , 2b): 

cr: T(2,26) ^ T(2,26) 

{x,y) {xPl mod 2 ,^). 

We have n 4 (cr((x,^))) = Wiix^y)) for all (x^y) G T(2,26). We say that u G I/^(T(2,26)) is antisym- 
metrie if uo a = —u and we denote by Aa the space of antisymmetric functions. 
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Let us note that for any function u on T(2, 2b) , —A{u o a) = (—Au) o a, so that Aa is stable under 
the action of — A. We write for the Friedrichs extension of the differential operator —A acting on 
C^(T(2, 26)) n Aa . The operator is self-adjoint, with domain iL^(T(2, 26)) fl Aa , and has compact 
resolvent. We denote by (A^’^)/c>i the sequence of its eigenvalues, which we will call antisymmetric, 
arranged in non-decreasing order and counted with multiplicities. 

Following m, we consider partitions of a specific type. For any positive integer £, we say that a 
-^-partition V = (T)^)i<^<^ of T(2, 26) is pair-symmetric if, for any i G {1,..., , cr{Di) = Dj with j i. 

Let us note that a pair-symmetric partition has an even number of domains, and that the nodal domains 
of an antisymmetric eigenfunction form a pair-symmetric partition. 

Lemma 2.5. Let V be a pair-symmetric and equispectral 2k-partition o/T(2,26) . We have 


xr<A2k{v). 


The proof is an application of the min-max characterization of eigenvalues, with test functions taken 
in the space Aa , as allowed by the definition of pair-symmetric partitions. We do not give the details 
here, and we refer instead the reader to m Proposition 6.3], where a similar result is discussed in the 
context of a double covering of the sphere. We now apply Lemma 2.5 to get a lower bound on the energy 
of a partition. 


Lemma 2.6. For 6 < \/\/k^ — 1, if V is a k-partition ofT{l,b) with no domain homeomorphic to a 
disk, then we have < A]^{'D) . 


Proof. A topological analysis of N{V) , using the hypothesis that no domain is homeomorphic to a 
disk, shows that has two distinct connected components for each i G {1, ,..., k} , which are 

exchanged by the map a . We refer the reader to HH for the details. Consequently, 114 ^{V) is a pair- 
symmetric and equispectral 2 /c-partition. According to Lemma [T5l we obtain 

< A2k{Ilf\V)) = Ak{V). 


A direct computation, using the facts that k is odd and that 6 < — 1, shows that . □ 


We can now complete the proof of Theorem 1.9 Let 6 < 6 ? and let P be a minimal /c-partition of 


T(l, 6 ). According to Lemma [2.4[ no domain of V is homeomorphic to a disk. From Lemmas 2.3 and 
| 2 . 6 [ we obtain 

A^(P^( 1 , 6 )) = k\^ < Ak{V) = (T(l, 6 )). 

This implies that P/c(l, 6 ) is minimal. 


2.3 Conjectures on the transition values 

Instead of the four-sheeted covering 1 X 4 : T(2,26) ^ T(l, 6 ), we now consider the two-sheeted covering 
112 • T(2,6) ^ T(1,6) , equipped with the map a : {x,y) i-G {x F 1 mod 2,y) . In the same manner 
as before, we can consider for this covering lifted partitions, antisymmetric functions, antisymmetric 
eigenvalues, and pair-symmetric partitions. In the rest of this section, these terms will be understood 
with respect to the covering 112 • T(2,6) ^ T(l, 6 ). Lemma 2.5 also holds in that case. We have the 
following conditional result. 

Proposition 2.7. IfV is a k-partition o/T(l, 6 ) such that Ilf^{V) is a 2k-partition o/T(2,6), then 


< A2k{Ii-\V)) = Ak{V). 


Proof Since we have assumed that is a 2 /c-partition, the pullback Ilf^{Di) of a domain of V 

has two connected components, and the map a exchanges them. This implies that Ilf^{'D) is a pair- 
symmetric partition, and the result follows from Lemma [tK] □ 


A direct computation shows that if k is odd and if 6 < 2/x/k‘^ — 1 , A^’^ = k‘^7r‘^ . If we were able to 
prove that a minimal /^-partition of T(l, 6 ) can be lifted to a 2 /c-partition of T( 2 , 6 ) when 6 < 2/x/k^ — 1, 
we would obtain bk > 2/\/k^ — 1. However, this is not obvious, even assuming that the boundary set of 
the partition does not contain any singular point (see m Section 5] for a discussion of this problem). 
We have on the other hand the following result. 
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Proposition 2.8. If k > 3 is odd, we have hk < 2/^/k‘^ — 1. 


To prove Proposition 2.8, we use the following result, whose proof is outlined in m Proposition 6.3] 
in the case of a double covering of the sphere. It consists in reproducing the arguments in the proof of 
na Theorem 1.17], while preserving the antisymmetry. 


Lemma 2.9. If V is a nodal 2k-partition assoeiated with an antisymmetrie eigenvalue, and if V has 
minimal energy among pair-symmetrie partitions, then A 2 /c(P) = 


Proof of Proposition 2.8. We have that P 2 /c( 2 , 6 ) is the partition lifted from P/c(l, 6 ), and is also the 
nodal partition of the antisymmetric eigenfunction (x, y) sin(/c 7 rx). Let us assume by contradiction 
that b > 2/\/^^ — 1 and is minimal. This wo uld i mply that V2k{2.,b) is minimal among pair- 

symmetric 2/c-partitions, and thus, according to Lemma 2.9, that = A2/c(T>2/c(2, b)) = , whereas 

a direct computation shows that the condition b > 2/y/k‘^ — 1 implies . 


□ 


Propositions 2T and 2 ^, and the numerical results of Section suggest the following conjecture. 
Conjecture 2.10. Let us denote b\ = 2/\/k2 — 1. If k > 3 is odd, we eonjeeture that bk = b\. 


Let us now try to analyse minimal partitions at the transition values. According to Proposition 


value IGtt^ . Since the eigenvalue IGtt^ has multiplicity 4, we obtain in this way minimal partitions 
which are not merely a translation of 1 ) 4 ( 1 , 1 / 2 ) . Figure shows an example whose boundary contains 
singular points. We conjecture that this partition is a starting point for the apparition of non-nodal 


1.8, 64 = 1/2. The minimal 4-partitions of T(l,l/2) are therefore nodal, associated with the eigen- 



Figure 1: A nodal 4-partition of T(l, 1/2) (associated with sin(47rx) + sin (47r^)). 


4-partitions of T(l,6) when b = 1/2 s, with 0 < 5 <C 1. More precisely, we conjecture that each 
singular point of order four splits into two singular points of order three (see Figures [7(d) [ 7(e)[ an dEl 
in Section 3.2.2 for numerical simulations). A similar deformation mechanism was already suggested by 
the numerical simulations in [2j Section 7] and [3l Sections 5 and 6], where the authors considered rect¬ 
angles and sectors, rather than tori, and where a singular point appeared on the boundary of the domain. 


In the case of an odd k , we are not able to give explicit examples of minimal /^-partitions which are not 
translations of P/c(l, 6). We can however cons truct candidates that would be minimal if Conjecture 2.10 
was true. For instance, for /c = 3, Conjecture 


2.10 


implies 63 = 1 /v^ and £ 3 (T( 1 , 1 /^/ 2 )) = Ott^ , which 
means that any 3-partition with energy Ott^ is minimal. We now look for antisymmetric eigenfunctions 
on T( 2 , l/\/2), associated with the eigenvalue Ott^ , which have 6 nodal domains. After projecting the 
corresponding nodal partition on T(l, l/\/2), we obtain a 3-partition with energy Ott^ . Figure]^ shows 
an example, in which the boundary contains singular points. In the same way. Figure shows how to 
obtain a 5-partition of T(l, 1/v^) by projection of a nodal 10-partition of T(2, 1 /V6 ), associated with 
the eigenvalue 257r^ . The former partition is minimal provided Conjecture 2 lOj is true. For k = 3 and 
k = b , the partitions obtained numerically, for b = 2/\/k‘^ — 1 -h £, seem close to these examples, with 
each singular point of order 4 splitting into a pair of singular points of order 3 (see Figures |5(d)| and 

m- 
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(a) A nodal 6-partition of T(2, l/\/2) (associated with (b) The 3-partition of 

cos(37ra:) — cos(7ra:) cos(2v^Try) + sin(7ra;) sin(2v^Try)). T(l,l/\/2) after projection. 

Figure 2: Construction of a 3-partition of T(l, l/\/2). 



(a) A nodal 10-partition of T(2, l/v^) (associated with (b) The 5-partition of 

cos(5Tra:) + sin(7ra:) sin(27r\/6y) — cos(7ra:) cos(27r\/6y))- T(l, l/\/6) after projection. 


Figure 3: Construction of a 5-partition of T(l, 1/v^). 


3 Numerical study 

3.1 Method 

For our numerical investigations, we adapt the method introduced by B. Bourdin, D. Bucur, and E. Oudet 
in [4]. In order to apply it, we approach the energy (1.1), interpreted as an infinity norm of the first 
eigenvalues, by a p-norm. 

Definition 3.1. For any 1 < p < oo and any k-partition V = (D^)i<^</c, we define 

A,,p(n)= ■ 


Then, we eonsider £k,p{M) = inf{A ; V G • 

In [4], the authors study the minimization of the sum of the first eigenvalues, which corresponds to 
the search for 2kp{M) in our notation. We extend the algorithm to cover the case p G [l,oo) and 
look for the minimal energy 2k^p{M) with 1 < p < oc large enough. This procedure is justified by the 
following result, proved in m- 

Proposition 3.2. The minimal energy £,k,p{M) is non-decreasing with respeet to p, and 


lim £fe p(M) = £fc(M). 

y^+oo 


To perform a numerical minimization, we mimic the method of [4]: we replace the minimal k- 
partition problem by a relaxed version, where we look for /c-tuples of functions (/i,..., fk) which satisfy 
Yli=i fi — ^ minimize a relaxed energy, depending on a small parameter 5 > 0. After performing 
a finite difference discretization of this problem, we work with a matrix 4> of size N x k. The integer 
N is the number of points in the finite difference grid, and the entry contains the (approximated) 
value of fi at the point indexed by /. We optimize a discretized version of the energy using a gradient 
descent algorithm. To ensure a better convergence, we start from a random initial data on a coarse finite 
difference grid, and we make progressive refinements 0 We refer the reader to [4] for details on all these 

^We thank Edouard Oudet for giving us detailed explanations on this point. 
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steps. In the end, we obtain a matrix ^ whose entries jire either 0 or 1 . We therefore have a discrete 
partition (i^i)i<i</c of the finite difference grid, where Di contains the points / such that = 1 . 

To give an approximation of £/c,p(T(l, 6)), from the result of the numerical optimization, we have 
two further steps which are not included in [4]: 

• construct a /^-partition V = of T(l, h) from (A)i<i</c; 

• compute the associated energy Ak{V). 

We want to construct a partition such that domains do not overlap and do not leave any part of 
T(l, 6 ) uncovered. Let us show how this can be achieved on an example. Figure [^represents a discrete 
3-partition of a 7 x 6 grid: the points of the grid are labeled by the domain to which they belong. We 
construct a boundary which separates points labeled by different integers. To represent this boundary, 
we could use a grid whose vertices are the midpoints of the initial grid. However, since we work on 
the torus T(l, 6 ), we can shift this grid by a half-step in the horizontal and vertical direction, and then 
deal with the initial grid. This gives us a strong partition of T(l, 6 ) (see Definition O- Then, we 
compute an approximation of its energy, without relaxation, using a finite difference approximation of 
the Dirichlet-Laplacian. 
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Figure 4: Partition obtained from the matrix 4>. 


Let us point out that this optimization algorithm is not always successful. All the other parameters 
being equal, it can converges rapidly to a good candidate for some initial data, whereas for others it 
terminates without reaching convergence. To overcome this problem, we have made several simulations, 
starting from different initial data, and compared the resulting energies. We present the best candidates 
obtained in this way. 


3.2 Results 

3.2.1 3-partitions of the torus T(l,6) 


We know from Proposition 


1.9 


2.10 


suggests that 63 equals 63 = l /\/2 


that 63 > 1/3. Conjecture 
0.707. If this is true, 1 ) 3 ( 1 , 6 ) should still be minimal for 1/3 < b < IIV 2 . To test this, we have 
implemented the method of Section 3.1 for b G {j/ 100 ; j = 34,..., 100 } . As was expected, the lowest 


energy in these cases is obtained for partitions of type 1 ) 3 ( 1 , b) as Figures [5(a)f[5(b)] show. 

Let us now study what happens when b is close to 63 . Figure shows our best result for different 
values of b . When b is greater than 0.71, the minimal partition seems to be a tiling of the torus by three 
isometric domains. These domains are roughly hexagonal, and close to the rectangles appearing in the 
partition of Figure [^ when b is close to 63 . For brevity, we will say in the following that a partition with 
hexagonal domains is hexagonal. 

For b close to 63 , the final result of the optimization algorithm appears to be very sensitive to the 
initial data. As a consequence, the partitions of Figures [5(d)| and |5(e)| were not actually obtained by 
starting from random initial data. Rather, we ran the algorithm, starting from a random matrix, in the 
case b = 0.81, where it produced an hexagonal partition similar to those of Figure We then used, as 
a starting point of the algorithm, the matrix obtained after two steps in the case b = 0.81. Of course, 
we compared our final results with those of other runs starting from random initial data, and found that 
they always had a lower energy. We used the same method for b close 64 when /c = 4, and for b close to 
65 and 1 when k = b. 
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(a) b = 0.64 



(b) b = 0.7 



(c) b = 0.71 


(d) b = 0.72 




(e) b = 0.73 (f) b = 0.8 (g) b = 0.9 

Figure 5: 3-partitions for some values of b. 


{h)b = l 


For a larger 6, up to 6 = 1, the best candidates produced by the algorithm are still hexagonal 
partitions, as seen on Figures [5(f)f 5(h)| For each 6, the energy of the best numerical candidate is an 
upper bound for £ 3 ( 7 ( 1 , 6)). This upper bound is plotted on Figure]^ as a function of 6, and compared 
with Qtt^ = As{'Ds{l,b)). We obtain a significant improvement for large b. The third upper bound, 
represented by a solid line, will be discussed in Section 



Figure 6: Upper bounds of £ 3 (T( 1 , b)) for b G {j/100 ; j = 30, ... , 100} . 


3.2.2 4-partitions of the torus T(l,6) 


We know from Proposition 1.8 that 64 = 1 / 2 . We are interested in the nature of minimal partitions for 
b close to 64 . Figure shows the best candidates obtained from the algorithm. 

As in the case /c = 3, we note the apparition of hexagonal partitions, shown on Figures [7(d)}|7(f) 


The hexagonal domains of the partition shown on Figure [7(d) | seem close to the square domains of the 
nodal 4-partition shown on Figure]^ This suggests that the partition of T(l,l/2) into four squares, 
shown on Figure is the starting point for the apparition of non-nodal 4-partitions of T(l, 6 ), when b 
becomes greater than 1 / 2 . 
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(a) b = 0.48 



(b) b = 0.49 



(c) b = 0.50 




(g) ^ = 1 


(d) b = 0.51 (e) b = 0.52 (f) b = 0.53 

Figure 7: 4-partition for b G {j/100; j = 48, ... , 53} and 6 = 1. 


For larger values of 6, up to 6 = 1, the minimal partitions are apparently still hexagonal. Figure 
|7(g)| shows for instance the best candidate for a minimal 4-partition of T(l, 1) . The energy of the best 
candidates gives us an upper bound for £ 4 ( 7 ( 1 , 6)), represented on Figure]^ as a function of b . 



Figure 8: Upper bounds of £ 4 ( 7 ( 1 , b)) for b G {j/100 ; j = 48, ... , 100} . 


3.2.3 5-partitions of the torus 7(1,6) 

We have conjectured that 65 equals 65 = l /\/6 ^ 0.408. The first images of Figure present the 
best candidates obtained numerically when 6 is close to 65 . They seem to support the conjecture. 
Furthermore, for 6 slightly larger than 65 , minimal partitions seem to be hexagonal, with domains close 
to the rectangles appearing in the partition on Figure 

For 6 between l /\/6 and 1, minimal partitions appear to be hexagonal. However, pairs of singular 
points in the boundary of these hexagonal partitions seem to merge when 6 approaches 1 , and for 6 = 1 , 
the best candidate produced by the algorithm is a partition of 7(1,1) into five equal squares. This 
process is shown on Figures 9(j) }9(1)[ We obtain again an upper bound for £ 5 ( 7 ( 1 , 6 )) as a function of 
6 , plotted on Figure pU) 
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I Biai^ 

(a) b = 0.40 (b) b = 0.41 (c) b = 0.42 


(d) b = 0.43 







(i) b = 0.9 


(j) b = 0.98 


(k) b = 0.99 


(1) ft = 1 


Figure 9: 5-partitions for some values of b. 



Figure 10; Upper bounds of £ 5 ( 7 ( 1 , b)) for b G {j/100; j = 18, ... , 100} . 


4 Examples of partitions 

4.1 Tilings of T(l, 6) 

The results of Section suggest that, at least for some values of k and b , the domains of minimal k- 
partitions of T(l, 6) are isometric polygons. In fact, except when k = b and 6=1, these polygons seem 
to be hexagons. On the other hand, according to Theorem |1. 2 [ any minimal partition satisfies the equal 
angle meeting property. 

This suggests the existence of partitions with a low energy, possibly minimal, that are tilings of the 
torus T(l, 6) by /c identical hexagons satisfying the equal angle meeting property. Let us note that this last 
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property is equivalent to the fact that the interior angle at each vertex of the hexagon is 27r/3. Finding 
these tilings is a purely geometrical problem. In this section, we prove Theorem |1.1Q| by constructing 
families of tilings, depending on k and 6, which seem close to the partitions obtained in Section We 
then compute the energy of these tilings, and compare it to the numerical results of Section 

4.2 Construction of the tilings 

4.2.1 Tilings of the plane 

To study tilings of the torus T(l, 6) by hexagons, it will be useful to consider tilings of the plane. As in 
Section [2^ we consider the natural projection map 

Hoc: ^ T(l,6) 

(x^y) 1 -^ {x mod 1,^ mod b ). 

We denote by (ei, 02 ) the canonical basis of M? , i.e. ei = (1, 0) and 02 = (0,1). 

Let us consider a strong and regular /^-partition V = {Di , ... , of the torus T(l, 6), such that 
all the Di^s are isometric to an hexagon that we denote by H . Let us note that, since V is strong, the 
area of H is b/k and, since V satisfies the equal angle meeting property, all the interior angles of H are 
27r/3. Let us then consider, for any i G {1, ... , /c} , the open set Il~^{Di) . It has an infinite number 
of connected components, each one being isometric to H . The family of all the connected components 
of all the sets , for i G {1,..., /c} , is a tiling of the plane M? by the hexagon H . This tiling is 

invariant under the translations associated with the vectors ei and be 2 . 

We can see that, conversely, the image by 11 of a tiling T of by an hexagon H is a regular k- 
partition of T(l,6) into domains isometric to H and satisfies the equal angle meeting property if the 
following conditions are verified by T. 

i. T is invariant under the translations associated with the vectors ei and be 2 . 

ii. The area of H is | . 

iii. All the interior angles of H are 27tI3 . 

We have therefore reformulated the original problem. We now look for the tilings of that satisfy 
properties i-iii, and, if possible, for an algorithm to construct those tilings. 


4.2.2 Change of basis 


Let T be an hexagonal tiling of 
useful to describe T. 


satisfying properties i-iii above. The following definition will be 


Definition 4.1. We say that a basis is T-adapted if its vectors connect the center of a tiling domain 
to the centers of two neighboring domains, with these two neighboring domains having a common edge 
(see Figure \Il(a)^ . 

Let us now denote by (ui,U 2 ) a T-adapted basis. 

Lemma 4.2. There exists a 2 x 2 matrix V = (vij) with integer coefficients such that 

vi,i , X 

Ui = —ei + —(6e2), 

Vl,2 . '^2,2/, X 

U2 = —ei + —(6e2), 


and det V = ±k . 

Proof. Since T is invariant under the translations associated with the vectors ei and be 2 , there exist 
integers 5iq, S 2 ,i, si ^2 and ^ 2,2 such that 


ei = Sipui + 52 ,iU2 ; 
be2 = Si,2Ui + S2,2U2 • 
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(a) 3-partition of a torus with (b) Area of the tiling domain, 

an adapted basis. 



(c) Construction of the tiling 
domain. 


Figure 11: Hexagonal tiling. 


With the notation S = (sij ), ui = tX2,i) ? ^2 = ('^^ 1 , 2 ,'^ 2 , 2 ), and U = (uij ), we have 


Thus detU det5' = b. But \detU\ is the area of the tiling domain H (see Figure p~l(b)[ ). Therefore 
detU = ±b/k and det5' = ±k. This implies that the matrix V = kS~^ has integer coefficients, and 
therefore has the desired properties. □ 

One can give a geometrical interpretation of the coefficients in the matrix V . Let us go back to the 
torus T(l,6) . We assume that the matrix V = (vij) is associated to the T-adapted basis (ui,U 2 ) (see 
Figure p"l(a) for an example). If we start from some hexagonal domain and translate it k times in the 
ui direction, it returns to its original position after turning uip times around the torus in the horizontal 
direction and U 2 ,i times in the vertical direction. Similarly, if we translate the domain k times in the U 2 
direction, it returns to its original position after turning vi ^2 times around the torus in the horizontal 
direction and ^ 2,2 times in the vertical direction. We can therefore say that the matrix V describes how 
the hexagonal tiling T wraps around the torus T(l, 6). 

The following result tells us at which condition we can solve the converse problem, that is to say find 
a tiling associated with a given basis. 

Lemma 4.3. Let (ui,U 2 ) be a basis , such that there exists a matrix V satisfying the properties 

of Lemma |^.4 Let Pq be some point in and let us note Pi = Pq + ui and P 2 = Pq + U 2 . The two 
following statements are equivalent. 

i. There exists an hexagonal tiling T such that (ui,U 2 ) is T-adapted. 

ii. There exists a point P in the interior of the triangle P 0 P 1 P 2 such that the segments PqP , PiP, and 
P 2 P meet with equal angles. 

Proof. Let us first consider the direct implication. We choose some domain of the tiling T and denote by 
Pq the vertex that connects the sides of this domain intersected by ui and U 2 (see Figure [Tl(c)| ). Then 
Pi and P 2 are vertices of the tiling, and there is another vertex P contained in the triangle P 0 P 1 P 2 . The 
point P has the desired property. 

Conversely, let us assume ii. After translating the three segments PqP, PiP, and P 2 P according to 
all the vectors of the lattice Zui + Zu 2 , we obtain the boundary of a tiling T which satisfies the three 
properties i-iii stated at the beginning of this subsection, and the basis (ui, U 2 ) is T-adapted. □ 

4.2.3 Reconstruction of the tiling domain 


We now try to determine at which condition Property ii. of Lemma [T3| is satisfied. We recall, without 
proof, a very classical geometrical result (see for instance [E]). 


Theorem 4.4. Let Pq , Pi, and P 2 be three non-eolinear points in (see Figure \i2(a)^ . One of the 
two following situations occurs. 
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(a) Equal angle property. (b) Construction of the Fermat point. 

Figure 12 : Fermat point of a triangle. 


i. If all three angles , ai, and 0^2 of the triangle P 0 P 1 P 2 smaller than ^ ^ there is a unique 
point P belonging to the interior of the triangle P 0 P 1 P 2 sueh that the segments PqP ; PiP and P 2 P 
meet with equal angles at P. The point P is ealled the Fermat point of the triangle P 0 P 1 P 2 . It is 
the point of minimum for the funetion Q QPq + QPi + QP 2 • 

ii. If ai > ^ for some i G {0,1,2}; then there is no point in the interior of P 0 P 1 P 2 at whieh the 
segments from the vertiees meet with equal angles. In that ease, the funetion Q QPq + QPi + QP 2 
reaehes its minimum at Pi. 


The following result gives an easy criterion for the existence of a Fermat point. 

Lemma 4.5. Let ui and U 2 be two non-zero veetors in M? . Let Pq be a point in M? . We set Pi = Pq + ui 
and P 2 = Pq + U 2 . The triangle P 0 P 1 P 2 has a Fermat point if and only if 


p G 




1 11 


2J 


or 


p G 


(^a) and p-\j 


1 — P‘^ 


< r < 


p-r-p 


(4.1) 

(4.2) 


with p = II II r = 

^ ||U1||||U2|| ||U2|| 


Proof The proof is a rather straightforward computation, and we merely indicate its steps. We express 
cosQ^i, for i G {0 , 1, 2} , as a function of p and r. Writing down the condition 


Vi G { 0 , 1 , 2 }, C0SQ!i G (- 5 ,1), 
we show that it is equivalent to the alternative (|4.1|) or (|4.2|). 


□ 


4.2.4 Algorithm 

The above results give us an algorithm to build a tiling of the torus T(l, b) by k hexagons: 
• choose a 2 X 2 matrix V with integer coefficients such that det V = zLk; 


• check whether the triangle generated by the vectors ui and U 2 (defined from V as in Lemma 4.2) 
has a Fermat point, using Lemma [45] 

• if the triangle has a Fermat point, compute its coordinates; 

• use the coordinates of the Fermat point to build the tiling domain. 
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Let us describe in more details how we perform the last two steps. We first recall a geometric 
construction of the Fermat point (see for instance [E]). 

Theorem 4.6. Let P 0 P 1 P 2 be a triangle in sueh that eaeh of the angles i G {0 , 1, 2} ^ smaller 
than ^ . Let us eonsider the three equilateral triangles lying outside 0/P0P1P2 and having one edge in 
eommon with it. For eaeh of these triangles, let us eonsider the line passing through the outer vertex and 
the vertex 0/P0P1P2 that does not belong to it (see Figure 12(b))^ . The three lines meet at the Fermat 
point P. 

Let us assume that we have performed the first two steps in the algorithm. We now have two vectors 
ui and U 2 such that the basis (ui,U 2 ) is adapted to some til ing T. We choose (arbitrarily) some point 
Po in . We can then build the points Pi and P2 of Figure 11(c) Using the construction of Theorem 
|4.6[ we can find the coordinates of the Fermat point P of the triangle P0P1P2 • The segments PqP , PP2 , 
and PiP then define three successive sides of the tiling domain, which is enough to construct the tiling 
domain itself. 


4.3 Examples 

We now look for examples of tilings with given matrices U, suggested by the numerical simulations in 
Section |3 
Theorem 


By applying Lemma [4.5[ we deduce the following result, which is a more precise version of 
TTOl 


Proposition 4.7. For /c G {3,4 , 5} , there exists a tiling of the torus T(l, b), with an associated matrix 
Vk ; if, and only if, b G (b^ , l] ; where 


• U3 = 


• U4 = 


• U5 = 




1 -1 
2 2 

1 -1 
2 3 


and b^ = ^; 


and 65 = 


H _ \/m-5\/3 
36 


If it exists, this tiling is a k-partition o/T(l, h), and we denote by H/c( 6 ) the corresponding tiling domain. 
We have 

£fe(T(l,6)) < min Ai(Hfe(6))), V6 e (5^, 1]. 

We use the finite element library Melina [18] to give an accurate upper bound of the first eigenvalue 
Ai(H/c(6)). This upper bound is represented in Figures [s] and [1Q| We observe that the upper bound 
by Ai(H/c(6)) is very close to the numerical simulations for b not too close to b\ (and not too close to 1 
when /c = 5). 


4.4 The case 6=1 

Let us consider the special case 6 = 1, in which the tiling domains can be described simply. The 
hexagonal tilings give a majoration of £/c(T(l, b)) , for /c G {3,4, 5} . Nevertheless, for /c = 5 , numerical 
computations show that a tiling of T(l, 1) by 5 squares, represented on Figure [T3(c) , has a lower energy. 

Proposition 4.8. We have the following upper bound on the minimal energy. 


i. £ 3 (T( 1 , 1)) < Ai(H 3 ( 1 )) 62.8389 where H 3 (l) is the hexagon shown in Figure 13(a) 




andH=^^-^^ 


with L = ^ , 


a. £ 4 (T( 1 , 1 )) < Ai(H 4 ( 1 )) 74.9467 where H 4 (l) is the hexagon shown in Figure 13(b), with L = ^ , 

^ ^ 4^ ’ = I “ 4^ ■ 

Hi. £ 5 (T( 1 , 1)) < Ai(Q) = IOtt^ 98.6960 where Q is a square of side ^ , and is the tiling domain of 
a b-partition ofT{l,l), as seen in Figure \l3(c) . 


Conjecture 4.9. The three inequalities in Proposition f.8 are actually equalities. 
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1 



(a) A hexagonal tiling domain for 
k = 3. 



(b) A hexagonal tiling domain 
for k = 4:. 


Figure 13: The case of T(l, 1). 



(c) A tiling by five squares. 


4.5 Comparison with the numerical results 


Let us focus for a moment on the case k = 3, for b close to 63 = l/\/2- Figure 14(a) 


shows that there 


exists 6 * > I/V 2 such that, for b G [^ 3 ,^*), Ai(H 3 ( 6 )) > Qtt^ . Therefore, for b G [^ 3 ,^*), the tiling of 
T(l,6) by three hexagons with straight edges that we have constructed is not minimal. The numerical 
method of Section generates better candidates. 




(a) k = 3 and b close to 63 . 


(b) k = 4 and b close to 64 . 




(c) k = 3 and b close to (d) k = 3 and b close to 1 . 


Figure 14: Zoom on the energy curves in Section 


This is consistent with the idea that there is some continuity of the minimal partitions with respect 
to b , and with the conjecture that the 3-partition shown on Figure [2(1^ is minimal. Indeed, if we try to 
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deform this latest partition by splitting each singular point of order four into a pair of singular points of 
order three, while keeping each domain close the original rectangle, the resulting partition cannot satisfy 
the equal angle meeting property if all the regular parts of the boundary remain straight lines. This 
suggests that the boundary of the partition should be curved in the neighborhood of the singular points 
in such a way that the we keep the equal angle meeting property. This seems to be the case for the 
partitions represented on Figure]^ This appears more clearly on Figure p~ 6 (a) , obtained by zooming on 
the neighborhood of two singular points. The same phenomenon occurs for 4-partitions of T(l, 6 ) when 

M(^ and | 16(b)[ an d fo r 5-par titions of T(l . 6 ) wh en b is close 


b is close to 64 = 1 / 2 , 
to 65 = 1 /or to 1 , 


as seen on Figures 
as seen on Figures 


the numerical method produces partitions with lower energy. 


14(c) 14(d), 16(c), and 15(d) As in Figure 14(a) 


we see that 



(a) b = 0.72, k = 3 (b) 6 = 0.51, k = 4 (c) b = 0.42, k = 5 W b = 0.99, k = 5 

Figure 15: Zoom on the boundary of some of the partition in Section 


For k G {3, 4, 5} and b G (^^,1] ? i^s denote by Tk{b) the hexagonal tiling of T(l,6) given by 
Proposition |4.7[ According to na Theorem 1.13], if Tk{b) is minimal, two neighboring domains of Tk{b) 
give a minimal 2-partition of the doubly hexagonal domain formed by their reunion. More explicitly, if 
we denote by 2H'^(6), j G {1, 2, 3} , the three doubly hexagonal domains we can extract from Tk{b) , we 
have £2(2H;^(5)) = Ak{Tk{b)) = Xi{Hk{b)). On the other hand, since we are considering 2-partitions, 
£ 2 ( 214 ;^(6)) = A2(2H;^(5)) . Therefore, if for some b G (6^, 1] , Tk{b) is minimal, we have 


Ai(Hfe( 6 )) = A 2 ( 2 H;^( 6 )) for j e {1, 2, 3} 

(this is a special case of [13 Proposition 8.3]). Figure [T^ gives, for k = 3,4, 5, the difference 
= Ai(Hfc( 6 )) - A 2 ( 2 H;(.( 6 )) for j G {1, 2, 3} and be{b'^,l], 


(4.3) 


computed thanks to the finite element library Melina [18]. It shows that in general. Equality (4.3) is 
not satisfied, and then the hexagonal tiling Tk{b) is not minimal. Let us note that when /c = 4 and 
b = \/3/2, the hexagonal tiling 71(5) is regular, and condition (4.3) is satisfied by symmetry. 



(a) k = 3 


(h) k = 4 


(c) k = 5 


Figure 16: b ha 5^(5), j G {1, 2, 3}, b G (5^, 1] for /c = 3,4, 5. 
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